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Abstract 

We present a functional renormalization group calculation for the two-dimensional (t, t')-Hubbard 
model at Van Hove filling. Using a channel decomposition we describe the momentum and fre- 
quency dependence of the vertex function in the normal phase. Compared to previous studies 
that neglect frequency dependences we find higher pseudo-critical scales and a smaller region of 
d-wave superconductivity. A large contribution to the effective interaction is given by a forward 
scattering process with finite frequency exchange. We test different frequency parameterizations 
and in a second calculation include the frequency dependence of the imaginary self-energy. We 
also generalize the channel decomposition to frequency-dependent fermion-boson vertex functions. 

PACS numbers: 71.10.Hf, 74.20.Mn, 75.10.Lp 



I. INTRODUCTION 



The phase diagrams of currently investigated quasi-two-dimensional correlated electron 
systems, such as the cuprates^, ruthenates^, and pnictides^i^, exhibit a variety of dynamically 
generated scales and phenomena caused by competing correlation effects. The fermionic 
functional renormalization group (RG) method has become a useful tool for the theoretical 
investigation of models for these materials, and more generally, low-dimensional correlated 
fermion models, both in and out of equilibrium (for a recent review, see Ref. |5|). 

The RG is a natural approach for understanding scale-dependent phenomena. Specific 
advantages of the method in its application to correlated fermion systems are that it does 
not require a priori assumptions about which correlations will eventually dominate, that it 
applies to a wide range of models, including multiband systems and models with general 
two- or even higher-body interactions, and that it allows to calculate the vertex functions 
of the associated quantum field theory, which are directly related to measurable quantities, 
as a function of momenta and (Matsubara) frequencies. Moreover, its application to weakly 
coupled Fermi systems is backed up by mathematical results implying full theoretical control, 
in the sense that for systems with weak, short-range interactions, convergent expansions of 
the RG map exist and can be used in the model class relevant for the above materials^"—. 
Practically, this means that the method is robust with respect to small changes and has 
good stability properties. 

In its most widely used form, the fermionic RG takes the form of an infinite hierarchy 
of integro-differential equations for the vertex functions^. Even in an approximation by 
truncation to a finite system of equations, it is difficult to take the physically important 
dependence of the vertex functions on momenta and frequencies accurately into account. In 
this paper, we report the results of a study for the two-dimensional Hubbard model in which 
we have taken into account the frequency and momentum dependence of the vertex functions, 
using the standard level-2 and Katanin truncation schemes described in Ref. Isl, as well as 
the vertex parametrization of Ref. Il3|. In the following we briefly review some important 
general arguments, issues, and results, first in the translation-invariant two-dimensional case 
and then in a more general context. 

Weak-coupling power counting indicates that in the RG fiow at low energy scales, the 
largest contributions to the fiow of the vertex functions come from the close vicinity of 
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the Fermi surface and small Matsubara frequencies. This suggests to start with a vertex 
function that depends only on the projection of momentum to the Fermi surface and of 
frequency to zero frequency as a first approximation. The vertex functions still depend on 
where the fermion momenta are on the Fermi surface. It is then natural to discretize this 
dependence, which converts the integro-differential equation into a system of many ODE. In 
the rigorous treatment of Ref. Q, the number of such points on Fermi surfaces with curvature 
is chosen to scale as A~^/^ as a function of the RG energy scale A. In practical applications 
to more general Fermi surface shapes, the number of patches is taken fixed, but as large as 
possible. This A^-patch method avoids any simplifying assumptions on the functional form 
of the projected vertex functions, and it has been very useful in understanding instabilities 
and drawing conclusions for phase diagrams^^""— . Its implementation at sufficiently large 
A^ is, however, numerically expensive for two-dimensional systems, since the number of 
components of a fermionic two-body interaction grows with the third power of A^. Using 
parallelization and present-day computing power, values of up to 144 have been attained, 
and the radial dependence has partly been taken into account^^. The parametrization of 
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approximates the four-point vertex by terms of the form of a boson exchange, which 
reduces the growth in in the resulting equations from cubic to linear. The price to pay 
is an additional error term, since representing the vertex function exactly in this way would 
require an infinite number of such terms. In the models we study here, however, it turns 
out that a finite, not very large number of them provides a good approximation in a large 
parameter ranged. 

The careful weak-coupling power counting argument is subtle, as it involves Taylor ex- 
pansions of the vertex functions and a detailed weighing of gains and losses in scale behavior 
in showing that projections can be taken. The argument is rigorous in the case of a single 
dominant instability and when the effective interaction is weak. However, the symmetry 
breaking at low temperatures manifests itself via a growth of the effective interactions for 
certain combinations of momenta as the RG scale parameter is decreased, so that the effec- 
tive vertex develops singularities corresponding to a slow decay of the effective interaction 
in position space. Although power counting improvements allow to continue the flow in the 
above-mentioned truncations^ beyond the region of small couplings, the truncation becomes 
unreliable already before a true singularity is reached. If one insisted on continuing the 
flow to arbitrarily low scales one would lose control over the errors. This so-called flow to 
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strong coupling is not a problem of the RG method as a whole, but can be mended when 
composite fields are treated appropriately. Then, the flow to strong coupling is found to in- 
dicate that correlations of composite flelds become long-range, leading to phase transitions. 
Since the dynamics of long-range correlations strongly depends on space-time dimension, 
the frequency dependence of the vertex functions becomes an essential issue in this regime. 
Besides allowing to check previous instability analyses, the results given here also provide 
vertex functions in a useful form for stopping the flow before the truncation fails, and for 
introducing bosonic order parameter flelds via a Hubbard-Stratonovich transformation. A 
good description of the frequency dependence is also important in the technical implemen- 
tation of RG flows for irreducible vertex functions, because their functional behavior at low 
scales is in general not simple, which raises the question how good simple parameterizations 
based on Taylor expansion around zero frequency really are. Secondly, the IPI equations 
never completely decouple the high-energy and high-frequency modes from the low-energy 
ones, so a quantitative description requires more than just the asymptotic low-scale infor- 
mation, but a good approximation in the entire frequency range. 

For the 2D lattice fermion systems, the frequency dependence was studied in various 
truncations. Besides being important for the analysis at low energy scales, tracking the 
frequency dependence of the vertex also becomes an issue when the frequency dependence 
of the fermionic self-energy is investigated. This is important in the study of quasi-particle 
weights and lifetimes. Inserting a frequency-independent four-point vertex in the standard 
equation for the self-energy flow in the IPI hierarchy gives a frequency- independent self- 
energy. In previous studies, the frequency-dependence of the self-energy was obtained in 
an approximation where bare propagators are used and the vertex flow is kept frequency 
independent, but the integrated flow equation for the vertex (which changes the one- loop 
frequency dependence) is inserted into the flow equation for the self-energy. In this way, 
the self-energy gets a frequency-dependence as in sunset diagrams, in which the momentum 
dependence is that of the flowing vertex. This approximation was used to calculate the flow 
of quasiparticle scattering rates^^, Z factors^ and to continue to real frequencies^"^ (Ref. 



2Juses an alternative formulation in the Wick ordered scheme). In a partially bosonized for- 
mulation, exchange propagators with a quadratic frequency dependence have been used^^i^i. 

In various further situations RG flows with frequency-dependent vertices have been stud- 
ied in some detail. In the case of the single-impurity Anderson model coupled to nonin- 
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teracting leads, integration over the degrees of freedom in the leads produces an effective 
local model in which all dependence is in the frequency variables. The dependence on the 
frequencies has been studied in great detail^^ and also at real frequencj*^. A representation 



3d . Very recently, the polaron problem in ultracold 



with boson fields was studied in Ref. 
gases was addressed^. 

This paper begins with a short recap of fermionic RG equations in Sec. [Tll In Sec. IIIII we 
describe our parameterization of the effective two-fermion vertex. The corresponding flow 
equations are derived in Sec. HVl In Sec. |V]we discuss different frequency parameterizations 
of effective boson exchange interactions. After numerical details are given in Sec. I VI I we 
present results at Van Hove filling in Sec. IVIIi 

We consider the extended {t, t')-Hubbard model on a 2D square lattice with unit spacing 

H[a\ a] = ^ ek4(k)a^(k) + U'^ n^,+n^,- , (1) 

where the kinetic term ek = — 2t(cos fc^ + cos ky) + At' cos kx cos ky originates from a tight 
binding approximation with nearest neighbor and next to nearest neighbor hopping. While 
t > sets the energy scale, t' G (0, 0.5)t measures the curvature of the Fermi surface 
{k G (— 7r,7r]^ : e(k) = ek — /i = 0}. We adjust the system to Van Hove filling, fi = —it', 
such that the free density of states is logarithmically divergent at the Fermi niveau. The 
interaction with [/ > is an onsite Coulomb repulsion of local densities n^^a = o)xa'^x,u- 



II. FERMIONIC RG 



We apply the fermionic Wilsonian RG based on integrating out fluctuations at scale A, see 
Ref. y for a review. A finite A > regularizes all graphs in the infrared, for example, cuts off 
momentum shells of thickness A around the Fermi surface. Our specific choice of introducing 
the scale A is different and explained below. We analyze the change of correlation functions 
when A is lowered in the RG flow within the level-2 truncation of the infinite hierarchy of 
flow equations^i2^. 
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The SU{2) and U{1) invariant two-particle vertex is parameterized as 

"1^a[^] = ^ j dpi . . .dp4 5{pi + P2 - P3 - P^) Va{pi,P2,P3) (2) 

X '^aiPl)^a'iP2)i^a'{P3)i^APA) , 
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where ijj and are anti-commuting Grassmann variables that replace the fermion operators 
of the Hamiltonian when writing the grand canonical partition function as a functional 
integral^^. The electron spin is denoted by cr and a' G {±} and the Pi = {p^, Pi) include both 
Matsubara frequency p^ and momentum pj. So the integrals / dpi are shorthand notation for 
J ^^^Ylpi thermodynamic limit with inverse temperature /3. One frequency and 

momentum is fixed because of energy and momentum conservation. The effective vertex 
function Va{pi,P2,P3) thus depends on three independent momenta and frequencies and the 
RG scale A. 

In our truncation the effective vertex function is determined by the initial condition 
Va{pi,P2,P3) = t/ for a — > oo and its flow equation^^ 

Va(pi, ■■■P3) = %Api, • • - ps) + ViXpu ■■■P3) + t;;{pi, ...P3), (3) 

where the dot denotes a derivative with respect to the scale A. The particle-particle contri- 
bution and the crossed and direct particle-hole contributions are respectively given by 

%piPi,---P3) = - j <^p[jk'^{p)G{pi+P2-p)\ Va{pi,P2,p)Va{pi+P2-P,P,P3) 

T^ipi,. . .P3) = - j dp [^G{p)G{p + ps - pi)] Va{pi,p + P3-Pi,P3) 

X Va{p,P2,P + P3-Pi) 
V^Pi^-'-Ps) = j ^p[jA^{p)G{p + P2-P3)\ '2Va{pi,p + P2-P3,p)Va{p,P2,P3) 

- Va{Pi,P + P2 -P3,Pl +P2 -P3) Va(p,P2,P3) 

-Vk{.Pi^P + P2-P3^p)Va{p,P2,P + P2-P3) , (4) 

where the scale derivative acts on both propagators but not on the vertex functions. Here 
G{j)) are fully dressed but regularized fermion propagators 

G{p) = -. TT^^TT' ^) 

tpo - e(p) - xa(p)Sa(p) 
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where Xa(p) is a multiplicative regulator function, which sets the scale A by suppressing 
frequency modes |po| ^ A or momentum shells with energy |e(p)| ^ A. For notational 
brevity, we have not indicated the A-dependence of G in the notation, but emphasize that 
both G and S (defined below) depend on A. For A — )■ the regulator function becomes 
unity so that G becomes the full propagator. The self-energy Sa(p) is independent of spin 
because of SU{2) invariance and can be calculated via 



Sa(p) = J dk S{k)[2VA{k,p,p) - VA{p,k,p)j, (6) 



where S{k) = G{k) — G{k)T,\G{k) is the single scale propagator, which is non-zero only 
where Xa(^) 7^ because 

= — 5; ~ '^^L • (7) 

(zfco - e(k) - XA{k)T.A{k)y 

Notice that we replaced single scale propagators by G in the flow equation for the effective 
vertex function. This is consistent in the level-2 truncation since S contains higher order 

Although on the right hand side of Eq. ([3]) only one-loop graphs appear, higher loop 
graphs are generated by solving the differential equation on iteration. The initial condition 
Va^oo{pi,P2,P3) = U is a, constant but the vertex function develops a singular momentum 
and frequency structure as A is decreased. A singularity in the flow signals the breakdown of 
the level-2 truncation and forces us to stop the RG flow at a stopping scale A*. By analyzing 
the momentum and frequency structure of the effective vertex at the scale A* we determine 
possible instabilities of the Landau Fermi liquid. Note that A* and the pseudo-critical scale, 
where a singularity appears in the vertex function, are only upper bounds of the critical 
scale for a possible phase transition. 

Already the bare one-loop bubbles $±(0 = j dp G{p)G{l ± diverge logarithmically 
for certain values of / in the limit A — )■ 0. At Van Hove filling the density of states di- 
verges logarithmically, which leads to an enhanced divergence of the particle-particle bubble 
$-(0) ~ (InA)^. Also the particle-hole bubble $+(/) diverges logarithmically at / = and 
/ = (0, 7r), with TT = (tt, tt), for t' /t G (0, |) at Van Hove filling. For finite but low enough A 
and t' 7^ 0, the regularized bubble $+(/) exhibits a maximum close to, but not exactly at, 
l = 7r. 
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III. CHANNEL DECOMPOSITION 



At not too low scales A ^ A* the vertex function is still regular and the potential 
singular structure of Eq. is given by transfer momenta pi + P2, Ps — Pi, and p2 — pj, 
that propagate through the fermion bubbles. In order to parameterize this frequency and 
momentum dependence we introduce three distinct channels^ 



^t{PuP4^P3-Pl) = T;,l{pi,P2,P3) (8) 
*k(Pi,P3,P2 - Pa) = -'2%i{PuP2,P3) + %"{pi,P2,Pl + P2 - P3), 

which we loosely call the superconducting, magnetic, and density channel, respectively. Each 
channel absorbs the graphs with one particular transfer momentum. A linear combination 
of these channels gives the flow of the vertex function, see Eqs. ([3]) and fl22l) below. The 
vertex function satisfies the symmetries 

Va{pi,P2,P3) = Va{ps,Pi+P2 -P3,Pi) 

= Va(p2,Pi,Pi +P2 -Ps) (9) 

= Va(Pi,P2,P3), 

that is, particle-hole symmetry, an anti-symmetrization property and time reversal sym- 
metry, respectively, and p = (— PO)P)- As a consequence of the first two symmetries the 
functions and $^ are symmetric in their first two momentum and frequency in- 

dices. This frequency and momentum dependence is assumed to be regular in each channel. 
The regularity is known to hold as long as the main momentum dependence is generated by 
a single channel in the fermion bubble. In case of competing instabilities at lower RG scales 
A this becomes a non-trivial assumption that, in the end, has to be checked. We will use 
this regularity to expand each channel in a sum of boson exchange interactions. 



In generalization of Ref. |13| we allow the fermion-boson vertex to depend on momentum, 
frequency, and the scale. For example, a general expansion of the superconducting channel 
is given by 
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The interpretation of the (m, n)-term in the sum is that two Cooper pairs corresponding to 
fermionic bihnears f dq ip_^_{q)r^{q, l)%lj_{l — q) couple in the form of a boson exchange, where 
the boson has propagator D^^{1). Here, m labels different gap symmetries (and possibly 
other details) of the pair. Note that in order to be a true boson propagator, D^{1) must 
satisfy positivity requirements, so that the Gaussian integral in a Hubbard- Stratonovich 
transformation in convergent. Because we stay within a fermionic formulation, we do not 
need to assume any positivity, so the ansatz is more general. In fact, not all terms come out 
positive in certain parameter regions (see below). For brevity, we shall still refer to them as 
boson propagators throughout. 

In this paper, we parameterize the fermion-boson vertex function 

Vi{q,l) = U\-c^)Tt{qoM (H) 

as a product of momentum-dependent form factors /„ that are orthogonal on [— 7r,7r]^ and 
boson-fermion vertex functions r^(go,^0;l) that depend on frequency and momentum. We 
choose a fixed orthonormal set {fn)n independent of scale. This would be no loss of generality 
if one retained infinitely many terms in the sum over m and n, but we shall keep only a small 
number of terms in this sum. This amounts to the assumption that the vertex function 
is regular. 

Apart from scale dependent coefficients D^^{1) we also allow each F^ to be scale depen- 
dent in order to capture part of the frequency dependence on go and q^. Due to particle-hole 
symmetry, i.e. the first line in Eq. ([9]), F„(go,^05l) = ^niJ-Q — for singlet supercon- 

ductivity, where /n(q) = /n(— q)- For simplicity, we expand in symmetric form factors only. 
A generalization to triplet superconductivity is straightforward since D^{1) consists of block 
matrices for singlet and triplet superconductivity. 

The scale-dependent remainder function 71^^ accounts for the fact that is not really a 
sum of finitely many boson exchange interactions. In the following we neglect this remainder. 
Ideally, the form factors /„ can be chosen such that 7^^^ contains only subleading terms. In 
practice, we want to deal with only a few expansion terms. How well the effective vertex 
function of the Hubbard model can be approximated by a few boson exchange interactions 
is analyzed in detail elsewhere^. 

Differentiating Eq. ( ITOl) with respect to the scale and integrating over two form factors 
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by using their orthogonality gives 

12„ j2_/ 



/ /md - q)/n(| - c{)^tc{<l, q', = Dt^mtiqo, k, l)r^{q'o, k, 1) (12) 



We use the normahzation T^{^-^,lQ,Vj = 1 such that the boson propagators are 



defined as in Ref. Il3|. For T = we set go = 9o = ^ (112) and obtain the fiow equation 

for the boson propagators in the superconducting channel 

DLii) = I (^^/™(i-q)/«(i-q')$^c(?,?',0|^^^^,^|. (13) 

By setting only Qq = ^-j but m = n in Eq. (fT2|) and using the fiow equation for we 
find the fiow equation for the boson-fermion vertex function 

rtiQO, lo, 1) = l^2§^2fni-2- q)/n(| " , |^,^| (14) 

|'?o=go="2" 

Inserting go = ^ gives r^(^,/o,l) = 0, which ensures that the normalization remains un- 
changed in the fiow. For T > the frequency ^ is not a fermion frequency in general. In 
this case we project to go = ±7rT and symmetrize over the sign. This corresponds, at T = 0, 
to a projection go = independent of Iq, which yields results very similar to the projection 
go = lo/2. 

Setting m = n is necessary since r„ does not depend on m and Eq. f|T2|) could otherwise 
not be satisfied. Later on we restrict to only two form factors describing s- and d-wave 
superconductivity. In this and other cases the boson propagators are nearly diagonaU^. We 
further use that the boson propagators have maxima in a neighborhood U around and tt 
in momentum space to choose only two values 1 = and 1 = vr for the dependence of F^ 
on the transfer momentum 1. We further restrict the frequency dependence to the relative 
frequency inside a Cooper pair. That is, we define 7^'^''*(go) := ^nilo^ 0, a) and write 

^tiqoJoA) = Yl 7f'^(| -go)i(/ G W(a)), (15) 

a=0,7r 

where W(a) = {p G (— 7r,7r]^ : \px + Py — a,x — ci,y\ < n} (periodically continued) are neighbor- 
hoods of a = 0, vr. The functions 7^'~^''* are symmetric due to time reversal symmetry (third 
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line in Eq. (|9])) and normalized to 7^°''*(0) = 1. Although we suppress this in the notation, 
the 7„ depend on the scale. Inserting into f|T^ and evaluating at /q = gives 

J2„ j2_/ 



*=o('(.9',«)l , „-7f'''to)*^cfe,?',<')l „, „ (16) 
with a = (0, a). 

In complete analogy we expand the magnetic and density channel in spin and density 
operators as well 

<{q,q'J)= E l(leW(a))5^Ml(/ + a)/„(^ + q)/„(^ + q') 

a=0,7r m,n 

X iTi% + 'ihTiq'o + 1 ) + Kiq, q', (i7) 

a=0,7r m,?i 

X 7:^'^(go + 'ihTiq'o + t) + g', /), (18) 

where one may choose different form factors than in the superconducting channel. As TZ^^ 
we also neglect 7^^ and TZ^ in the following. For the boson propagators we obtain 

^1(0= /7f^^/md + q)/n(| + q')<^>^:;(g,g',0l , k (19) 

kL{i) = J (^1)2 ^2n^2 U\ + q)^(l + qO^;^(g, q\ 

and likewise for the frequency dependent part of the fermion-boson vertex 
7r(go) = / /„(q)/„(q') [$;!:(?, g', a) |^,^^-7r(go)$;::(g,g', a) 



90=<?o=0. 

(20) 



IV. DETAILED SET-UP AND FLOW EQUATIONS 

We choose a minimal set of form factors that describe the leading instabilities of the 
{t, t')-Hubbard model at Van Hove filling, and at not too large t' for all fillings. In the 
superconducting channel f\ (p) = 1 describes s-wave superconductivity with a peak at 1 = 
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and s-wave alternating pairing with a peak near 1 = vr. Both couphngs are subleading 
but are essential at higher scales, and ultimately suppress the pseudo-critical scale to lower 
values. A second form factor /2(p) = cosp^: ~ cos Py in the superconducting channel is 
used to describe d-wave superconductivity. In the two particle-hole channels we only use 
/i(p) = 1- Then a singularity of M^(/) at / = signals ferromagnetism and a singularity 
near (0, vr) (incommensurate) antiferromagnetism. 

In principle -D^2(0 couples s- and d-wave superconductivity. The coupling, however. 



is small and exactly zero for |/^.| 



\ly\ at Iq 



0, which includes the peaks at and vr. 
We therefore neglect this coupling and consider only diagonal boson propagators. In order 
to parameterize the momentum dependence of boson propagators B = D,M or K we 
distinguish their values close to and n 



(21) 



a=0,7r 



where from now on we drop the explicit notation of the scale dependence. 

The flow equations for the boson propagators, fermion-boson vertices, and the fermionic 
self-energy can now be derived by inserting the channel decomposition of the effective vertex 
function 



Va{Pi,P2,P3) = U - $^c(Pl.P3,Pl +P3) + $m(P1'P4,P3 -Pl) 

+ ^*m(Pi>P3,P2 - Pa) - ^*k(Pi>P3>P2 -ps) 



(22) 



in the equations of the last section. We use that the boson propagators are symmetric in 
their dependence on each momentum and ly separately. In the superconducting channel 
we arrive at 

d 



dp 



dA 



G'(p,|-po)G'(l + a 



P,|+Po) 



fn 



l+a 

^ 2 



P) 



(23) 



^:(07f'>o)-f/5n=l-«f(PO,/o) 



1 2 



with the following contributions to the frequency dependence from box and vertex diagrams 



(Po,/o) 



W(0) 



d^k 

(27r)2 



A„(k) 3M°(k,po)7°(^)7°(^; 



(24) 



-l)"-i3Mf(k,po)7^(^)7^(^: 



-(-1 



3AT(k,po)7^(^)7^(^^ 
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with Ai(k) = 1 and A2{k) 



cos kx +COS ky 



for s- and cZ-wave superconductivity. Notice that 



the projected momentum dependence is the same for direct, box, and vertex diagrams and 
is just given by a form factor. The function A2 is obtained by using simple trigonometric 
identities and the symmetric momentum dependence of the boson propagators mentioned 
above. 

The flow of the frequency dependent part of the boson-fermion vertex function is given 

by 



7r'^(go) = 



dp 



^Gip)Gia - p)] /„(§ - p)2 iDmin'^Po) - U5n=i - «f (po, 0)' 



dA 

X 



'U5n=, (7r'^(go) - 1) + 7r'^(go)cif (Po, 0) - if (po, go)] (25) 
with frequency dependent (unprojected) box and vertex contributions 



if bo, <lo) = \j 7f^^n(k) [3M?(k,po - qom^f 
^ Ju{0) [^'^) L 



(26) 



-l)"-i3Mf(k,po-go)7r^"^' 



-K?(k,po-go)7K 



0(Po+go\2 



- (-1 



Notice that A^{po,0) = a^{po,0) represents the frequency projection. 
Similar flow equations are obtained for the magnetic propagator 



Mm 



dp ^G{p,po-^i)G{p + \ + a, 



Po + I^ 



(27) 



U + Mf (/)7^,(po) + ar(po, lo) + a"(p, I + a) 



with 



a=0,7r 



a=0,?T 



(27r)2 L 'Mv 2 ^ iM\ 2 



Mf(k,po) -7K(^)7K(^)i^r(k,Po) 



2 7s\i(^)7sai(^)^r(k,Po)] 



(28) 



d^k 



(2 



TT 



2 7sc,2V 2 'SC,2V 2 



)7 



D2^(k,Po)/2(p-f)/2(p + l-f) 



Using trigonometric identities the momentum dependence of a^(p, /) can be made explicit. 
The flow equation of fermion-boson vertex function reads 



7m (90) 



Mf(0) 



dp 



—G{p)G{p + a) 



U + Mf (0)7^(po) + <ipo, 0) + a^{p, a) (29) 



X 



U[l- 7m('?o) + A^po, go) + A^ip, a, go) - 7m(?o) a^Po, 0) + a^(p, a) 
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with 



<(Po,go) = ^ E / [7S(2^)2Mf(k,po-go) -75(^)'^f(k,Po-go) 



A^{p,a,qo) 



-2 7sc,i(^)'^r(k,Po + go 



(30) 



a=0,7r 



(27r)2 '^'^•^^ 2 

Similar flow equations are obtained in the density channel 

k-,(i) = - [dp 



7sa.(^)'A"(k,Po + go)/2(p - |)/2(P + a - I). 



d 

dA' 



T^G'(p, po - ^^(p + 1 + a, + ^) 



(31) 



X 



K^ilhuiPo) -U- a^ipo, lo) - a^ip, I + a) 



with 



af(po,/o) 



d^k 



2 ^ J (27r)2 



[37M(^)7M(^)Mf (k,po) + 75(^)7K(^)^^(k,Po) 



2 7sc,i(2^)7sc,.(2^)^r(k,Po)]. 



(32) 



The boson-fermion flow is given by 



75(go) 



Km 



dp 



—G{p)Gip + a) 



iqiOh^iPo) -U- afipo, 0) - afip, a) (33) 



U{1- 7k(9o) + A^ipo, go) + ^"(P, a, go) - 7^(go) af (Po, 0) + < (p, a) 



with 



a=0,7r ' 



d^k 

(2^ 



[37S(2^)'Mf (k,po - go) + 7S(^)'^f (k,Po - go) 



2 'y^^ ^^2!lzJo ^2 na 



z^f(k,po + go)] • 



(34) 



Suppose that Xa 7^ for all finite A, that is, arbitrarily large scales have to be taken into 
account. In principle, the initial condition could be stated for A — j- oo as the bare action 
coming from the Hubbard Hamiltonian. We find it more convenient to choose a relatively 
large scale Aq ^ U and treat all scales A > Aq by perturbation theory. This is possible 
because Aq is an infrared cutoff, such that the small parameter is effectively given by 
The effective vertex in second order perturbation theory reads 



^Ao(Pl,P2,P3) = U- f/2<l>^^°(pi +P2) - f/'$^°(p3 - Pi) + 0{^J' 



(35) 



14 



with contributions from the particle-particle and the particle-hole channel. We read off the 
initial condition for s-wave superconducting boson propagator 

Dlo{l) = ^^'^{l + a) (36) 

and choose for the magnetic and density boson propagators 

Ml,{l) = Kl,{l) = <l>lo{l + a). (37) 

Since we treat both particle-hole channels in the same approximation, a different assignment 
of the initial condition among the magnetic and density channel would not influence the 
outcome of the flow. In all these s-wave channels the initial boson-fermion vertex function 
is equal to one, 

7r'^M=7MM=7K(go) = l (38) 

for a = or ^. 

Initially, there is no coupling for rf-wave superconductivity, so D|o(0 — 0) where for the 
remainder of this section we only consider a = in the d-wave channel and suppress the 
explicit notation. Therefore the fermion-boson vertex function is undefined in this channel. 
In order to circumvent ^ in Eq. fl2^ we propose two solutions 

1. We arbitrarily set 7|'^(5'o) = 1 and compute the flow with constant 7!^^ from Aq down 
until .02(0) reaches a minimal value D^^^. This defines a scale Ai, where we switch on 
the flow of 7!*^. The flow rapidly adjusts from a constant 72°(q'o) = 1 to a decaying 
function. We checked that the result does not depend on reasonably chosen minimal 
values G [10-^10-3]. 

2. In order for ^110^02(0)^0 ^2'^ {%) to exist it is necessary that 

hm D2(0)7f (go) = 0. (39) 

-D2(0)->-0 

This condition is satisfied if 

sc. I dp [A.G{p)G{-p)] /2(p)^af (po,0)if (po,go) 
^ Idp[£G{p)G{-p)]MpyaT{Po,Oy ' 

Taking this as the initial condition ensures that 71^^ (go) has the same frequency decay 
as the vertex and box diagrams (but normalized) already at the beginning of the flow. 
This in principle allows to start the full flow at Aq with D2{1) = and 7|'^(go) as 
above. 
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For the time being we have apphed the first variant for simphcity. Since the effect of 
non-constant fermion-boson vertex functions turns out to be not very large we resign from 
implementing the second variant, which is fundamentally more sound. 

The flow equations and initial conditions can be evaluated for different regulator functions 
Xa(p)- Physically, the result should not depend on the choice of regularization. However, 
because of the level-2 truncation and the requirement to stop the flow once singularities 
appear there is a cutoff dependence of quantitative results. In order to check our results 
qualitatively we use two very different regularization schemes. At temperature zero we use 
a soft frequency regularization. Here the RG scale A is denoted by fl, and 

In a second scheme we use temperature A = T > as a regularization as introduced in Ref. 
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Here the Grassmann fields \& are rescaled such that the quartic part of the microscopic 
action does not explicitly depend on T. This results in a regularized free propagator that 
depends on T through the discrete frequency variables and an overall factor ^/T. 

To complete the set-up of our study, we discuss our determination of the 'transition' scale 
f2* in the fi-scheme, and the 'transition' temperature T^, in the T-scheme. As in previous 
studies, we never run a flow up to, or even very close to the point Qq (resp. T^) where any 
part of the vertex function is divergent, because the truncation we use breaks down before 
that point. Instead we define a stopping condition, namely that the largest value of one of 
the exchange propagators reaches Vmax = 20t, which corresponds to about 2.5 times the free 
bandwidth. While the choice of Knax is arbitrary to a certain extent, we have verified that 
varying Vmax > 16t does not change our results qualitatively. 



V. TRANSFER FREQUENCY DEPENDENCE 
A. Lorentzians as Vertex Functions 



In this section we motivate an explicit frequency parameterization of the boson propaga- 
tors. The goal is to capture the leading frequency behavior for small scales near the critical 
scale. 

Consider the case where one channel is clearly dominant, for example, a curved and 
regular Fermi surface, for which Umklapp scattering can be neglected. There, the onset 
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of superconductivity determines the stopping scale. Neglecting the marginal particle-hole 
graphs corresponds to neglecting vertex and box diagrams in the flow equation for the 
superconducting boson propagator, which can then be solved to give the RPA result 

" 1 - ' ^^^^ 

where = J dp G{p)G{l — p)fn{\ — p)^ is the particle-particle bubble and -D°„(/) 

is the initial interaction. For notational simplicity we restrict to only one form factor. 
The following analysis generalizes to more form factors since the boson propagators are 
approximately diagonal. 

Suppose that D^^{1) = u is a positive constant, that is, an attractive initial interaction. 
For momentum 1 = 0, small scales fi, and small frequencies /q the particle-particle bubble 
at Van Hove filling is roughly given by 

Re $r (/o, 0) ~ +C In^ {^^ + Zll) . 

where C and Z are positive constants and we have neglected imaginary parts of the bubble. 
That is, the bubble diverges for f2 = as In^ |/o|- The leading singular structure of D{1), 
however, does not contain a logarithmic frequency term. Instead a singularity at a scale 
VLq builds up. It is characterized by a mass term = 1 — uClri'Vl?' that vanishes at this 
critical scale VLq. Just above the critical scale the denominator of Eq. (l42l) at 1 = can be 
written as 

1 - uCln'(fi2 + Zll) = nil-uC\n(^ + Z ^ (^ln(fi2 + Z/^) +lnfi2j 

= ml-AuCZ\nn-(^^ +C((/oM)^), > l^c > 0. 



A similar argument holds for the particle-hole channels, where the particle-hole bubble 
diverges with Infi for / = at van Hove filling. Due to the fi-regularization there is no 
damping term proportional to ^ for small /qj and at > VLq. 

The imaginary part of the particle-particle bubble vanishes at zero frequency. However, 
it exhibits singular behavior in its first frequency derivative if the density of states is not 
symmetric. At Van Hove filling and for curved Fermi surfaces, < t' < t/2, 

(9iolm$i^)(/o = 0,0) = -Ciln2l] + C(lnl]), Ci > 0. (43) 
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This kind of small frequency behavior can be captured by parameterization of the fre- 
quency dependence of the boson propagators in each channel with "Lorentzians" 

Bmrn{lo,\) = [mli\)+tanm + bUl)lir' (44) 

with momentum dependent parameters mQ,an,bQ G M. A flow equation for these param- 
eters is obtained by considering the RG equation for Bmmi^o^ 1) and its first two frequency 
derivatives at zero frequency /q = 0. This corresponds to a second order derivative expan- 
sion with multiple bosonic fields. Parameters are thus determined from the small frequency 
behavior only. In addition to small frequencies also the Iq^ decay of s-wave channels for 
large frequencies Iq is described correctly but with a wrong prefactor. The d-wave super- 
conducting boson propagator is solely generated by other channels that decay in frequency 
themselves. It therefore has an Iq"^ decay at finite 1] > 0. 

In RPA, the flow equations are decoupled such that the leading instability is determined 
by the evaluation of direct graphs in the channel decomposition only. Here only the singular 
value of boson propagators at frequency zero determines the instability. In contrast, the 
fermionic RG couples different channels via box and vertex diagrams. The value of these di- 
agrams is computed by a convolution in loop and vertex frequencies. Therefore, the coupling 
of different channels may be estimated wrongly when intermediate to large frequencies are 
described inadequately. In the next section we compare Lorentz curves to the full frequency 
dependence of boson propagators calculated in the RG flow. 

B. Comparison to Actual Dependence 

The Lorentz parameterization of the transfer frequency has been motivated above for the 
small frequency regime. There is no obvious reason why it should extend to an accurate 
description of exchange propagators for arbitrary frequency values. This is why, in a second 
calculation, we discretize the transfer frequency dependence and then trace it fully during 
the flow. This allows us to check the Lorentz ansatz. 

In anticipation of the full discussion of the RG flow with frequency dependent vertices we 
show in Fig. [1] several approximations to the discretized frequency data for a typical example, 
the antiferromagnetic exchange propagator at the stopping scale Q^,. By construction, a 
single Lorentz curve determined from numerical extraction of Taylor coefficients reproduces 
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FIG. 1: (Color online) Frequency dependence of the magnetic propagator Iq i— t- Mii(/o,7r) for 
parameter values U = 3t, t' = 0.2t at the stopping scale ri* = 0.089 t. Plus marks indicate 
discrete frequency data from the RG calculation. A simple Lorentz parameterization from the 
small frequency expansion (dashed line) accurately describes the small frequency behavior. The 
least squares fit to the sum of two Lorentz distributions (solid line) also captures well the frequency 
dependence away from the small frequency region, it is composed of the two dotted curves. 

the discrete reference data well for small frequencies (dashed line in Fig. [T]). However, this is 
not true for frequencies /q > At Iq ~ 2^2* the relative error is about 15% and it becomes 
very large at higher frequencies. We find a similar situation in the temperature flow setup, 
this scheme is shortly discussed in Sec. IVII CI Fitting a single Lorentz curve to the discrete 
data with a least square condition brings no improvement. 

Compared to other RG schemes, in the IPI scheme contributions to the flow at scale A 
are not restricted to low frequency or low energy processes with |e| , < A. Moreover, 
the f2-scheme provides just a mild regularization of the free propagator such that its scale 
derivative as well as the corresponding single scale propagator do not have compact support 
in frequency space. Hence, proper frequency parameterization away from the small frequency 
region possibly is important. 

Such a parameterization can be accomplished by a sum of two Lorentz curves (solid line 
in Fig. [T]). One of these curves can be interpreted as a small frequency process and the 
other one as a large frequency process. Parameters in this ansatz have been determined by 
a weighted least-squares fit of the discrete data, with a ten times larger weighting factor for 
frequencies with \uj\ < Q. 

In Sec. IVIII we compute the flow with a single Lorentz curve and compare to the flow 
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(tt, tt) 



(0,0) (0,7r) 

FIG. 2: Illustration of momentum discretization for exchange propagators in 12 segments about 
k = (0,0) and k = (vr, vr). Small circles mark the representative momentum associated with each 
segment. In calculations we use an analogous discretization in 32 segments. 

with discretized frequencies. Surprisingly, the simple Lorentz curve ansatz captures well the 
flow of the leading couplings of the interaction vertex in most parameter regions. However, 
it does not detect a new type of scattering singularity, detailed below, that we encounter in 
the setup with frequency discretization. Also it overestimates the flow of the Z factor at 
momentum (0, vr) for all scales and produces scale derivatives Z with a relative error of 
up to 40%, see section IVTlDl 

On the other hand, we find that quantitatively accurate results can be obtained with the 
frequency parameterization by a sum of two Lorentz curves. Although its effectiveness for 
the extraction of information about arbitrary observables is not tested, this parameterization 
reproduces very well the flow of the interaction vertex as well as the flow of the Z factor. 

This functional form may turn out to be a promising generalization of bosonic field 
theories for future applications in condensed matter theory. A simple strategy for determing 
its parameters, independent of a prior frequency discretization procedure, needs yet to be 
established. 

VI. NUMERICAL IMPLEMENTATION 
A. Transfer Momentum Parameterization 

The momentum dependence of boson propagators is discretized using step functions. By 
reflection symmetries about the coordinate axes and the Brillouin zone diagonal one can 
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restrict to discretization of one eighth of the Brillouin zone, see Fig. [2l The discretization in 
radial segments about k = (0, 0) and (vr, vr) permits a detailed momentum resolution in the 
vicinity of these points. In case of incommensurate antiferromagnetism the maximum of the 
magnetic boson propagator can move far away from (11,71) in the (1, 0)-direction. Thus we 
place several representative momenta along the coordinate axes. In the calculations below 
we use a discretization in 32 segments per one eighth of the Brillouin zone for each exchange 
propagator. 

B. Frequency Discretization 

Within each momentum segment, the frequency dependence of exchange propagators is 
discretized. By symmetry, consideration of transfer frequencies Zq > is sufficient. We use 
a logarithmic grid in the frequency range 10^^ < lo/t < 3 ■ 10^ and include the frequency 
value lo = 0. 

We have checked that our numerical results are not affected by the concrete implemen- 
tation of frequency dependences. For this purpose we have (a) used a frequency grid of 20 
points, which still allows application of adaptive numerical integrators to the loop frequency 
integral on the rhs of flow equations, and (b) used a frequency grid of higher resolution, with 
evaluation of frequency integrals as a discrete sum. 

The fermion-boson vertex functions 7(^0) ^-re symmetric around go = as well. We 
choose a frequency grid of 34 base points for positive frequencies. For the evaluation of the 
flow equations the boson-fermion vertex functions have to be evaluated at qq, "^""^^ , and 
30^0 ^ where po is the loop frequency and go the external frequency. The latter two values 
are obtained by linear interpolation. 

C. Loop Integrals 

The evaluation of loop frequency integrals in IPI flow equations by contour techniques 
is still possible in the presence of frequency dependent vertex functions but becomes more 
involved. In calculations where we use a momentum independent cutoff and neglect the self- 
energy the evaluation of loop integrals can be simplified by first integrating over the loop 
momentum and then over the loop frequency. The crucial point is here that the momentum 
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integrals are independent of the RG scale and have to be evaluated only once. More precisely, 
the free propagator with Q cutoff Cn{p) = Xn{po) {ipo — ^p)~^ allows the factorization 

^(Cn{p)Cn{k)') = Co{p)Co{k) ±(^xn{po)xn{ko)) (45) 

in a momentum independent and a scale independent term. 

The flow equations for exchange propagators have the general form 

B^mil) = / ^ ^(^^(^° - 2)X^(^'0 + I)) E"?i(^0,P0)«g(/0,P0) 4n,^Jl,P0), 

4,,^,,(^>Po) = C'o(±(p- |))Co(p+ I) ^,,(l,p)^,,(I,p). (46) 

Here the dependence of the projected interaction vertex on loop and external momenta 
is extracted analytically using trigonometric identities^. This produces two sums of fre- 
quency dependent functions and momentum dependent functions ipj in (1461) . where ipj 
is scale independent. Therefore, the momentum integrals I^.^^^.^{1,pq) are scale indepen- 
dent and can be calculated before starting the flow. This needs to be done for all discrete 
frequency- momenta / = (/q^I) where exchange propagators are calculated and in addition 
for all combinations {-^^ipj-^^ipj^) that occur. 

We discretize the momentum integral expression in the frequency variable po as follows: 
By symmetry, only Pq.Iq > is needed. ^, i^,Po) is regular in pq except for possible 
singular behavior at po = /o/2, which drives the flow. The asymptotics of the momentum 
integral for zero external momentum 1 = and loop frequencies po close to Iq/2 can be 
explicitly calculated. At Van Hove filling and with ipj-^{\,p) = ipj^iljp) = 1 we find 

/A((/o,0),po)= ''''''' (47) 
[-C2Sgn(po-/o/2)lnbo-/o/2| + 0((po- V2)°) : /o > 0, 

with positive constants Ci, €2- Typically, I^, ^. {l,Po) is a monotonic function in po at both 
sides of the singularity point po = lo/2. 

To capture the structure of such momentum integrals we discretize the frequency variable 
Po on a grid that depends on Iq. We use logarithmic spacing and about 100 data points at 
both sides of the point po = Iq/2. The minimal distance to this point in the grid is adjusted 
depending on the expected scale usually f2*/10 is sufficient. This ensures that during 
the flow the full structure of integrands in the RG equations is properly taken into account 
at all scales. 
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VII. RESULTS AT VAN HOVE FILLING 



A. Results in the Q-Scheme 

In this section we discuss several frequency parameterizations of the interaction vertex, 
going from simple setups to more detailed parameterizations, and compare the resulting flows 
for systems at Van Hove filling and temperature zero. Here, the self-energy and imaginary 
parts of the vertex function are neglected at first, both will be discussed in the next sections. 

In the left graph of Fig. |3] the stopping scale is plotted over next to nearest neigh- 
bor hopping t' at U = 3t. The lowest line corresponds to frequency independent vertex 
functions. This approximation is motivated by weak coupling power counting and has been 
used in early RG applications to the 2D Hubbard model^"— . In agreement with Ref. |l9 



we find regions of antiferromagnetism (AFM), d-wave superconductivity (d-SC), and fer- 
romagnetism (FM) at Van Hove filling. That is, for low t' the dominant coupling is the 
boson propagator in the magnetic channel with form factor /i(p) = 1 at (commensurate 
AFM) or near (incommensurate AFM) momentum transfer (vr, vr). For intermediate hopping 
ratios t'/t the fermionic interaction has a sharp peak in the superconducting channel with 
(i-wave form factor /2(p) = cosp^; — cos Py at zero transfer momentum. For ^ < t'/t < ^ 
the ferromagnetic instability is signaled by a divergence of the magnetic boson propagator 
with form factor /i(p) = 1 at zero transfer momentum. In the crossover region between 
superconductivity and ferromagnetism the stopping scale is strongly suppressed compared 
to t' = 0. 

Our simplest frequency parameterization considers frequency independent fermion-boson 
vertex functions, r^(go,^0!l) = 1? and describes the frequency dependence of the boson 
propagators with one Lorentz curve. The stopping scale obtained in this approximation is 
the upper most curve (ii) plotted in Fig. [3l There are two major differences to a frequency 
independent approximation, line (i). First, the stopping scale is much higher, especially for 
intermediate next to nearest neighbor hopping. Secondly, there is no region of dominant 
d-wave superconductivity anymore. Instead, incommensurate antiferromagnetic couplings 
with transfer momenta quite far from (tt, vr) are dominant for intermediate t'. 

In our understanding both effects are closely related. In general, the static approximation 
overestimates the vertex function. Consideration of a frequency dependent vertex function 
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FIG. 3: (Color online) Stopping scales for systems at Van Hove filling. Left: Comparison of 
different frequency parameterizations: (i) Frequency independent vertex, (ii) Vertex with simple 
Lorentz parameterization of the transfer frequency, (iii) Vertex with discretized transfer frequency, 
(iv, symbols) Vertex with discretized transfer and non-transfer frequencies. Right: Setup (iii) 
for interaction parameters U /t = 3, 2.5 and 2 (top-down). Most dominant ordering tendencies: 
commensurate AFM (dashes, filled square), incommensurate AFM (dotted, open square), d-SC 
(solid, circle), scattering instability (dash-dot, triangle), FM (short dashes, diamond). 

leads to an additional decay of the frequency integrand in box and vertex diagrams, unless 
they are generated directly by U . The latter is the case for all s-wave channels with form 
factor /i(p) = 1 but not for the (i-wave superconducting channel. Therefore, box and 
vertex diagrams that generate an attractive (i-wave interaction contribute less after the loop 
frequency integration. More generally, the coupling of different channels is reduced due to 
the frequency dependence of box and vertex diagrams since they are not evaluated at their 
maximal value only. The effective width of boson propagators approaching a singularity 
tends to zero, but is also reduced for all other boson propagators at low scales, see Fig. HI This 
leads to a reduction of screening and mutual coupling between the channels and consequently 
to a higher pseudo-critical scale. 

We stress that d-wave superconductivity is still generated by the RG flow, but to a much 
lesser extent. Due to the higher stopping scale this generating process has not enough 
RG time to become a leading instability. If we change our definition of to allow lower 
scales then rf-wave superconductivity becomes dominant in a small parameter region of t' 
eventually. For example, this can be achieved by not taking the maximal value of the 
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dominant boson propagator at frequency zero as the stop condition but rather a frequency 
mean value over a not too small neighborhood around zero. (This would not change f2* in 
case of a frequency independent calculation.) However, the so defined stopping scale is still 
much larger than the stopping scale obtained with frequency independent vertices and the 
region of rf-wave superconductivity is substantially smaller. 

For lower initial interaction values U the effects of frequency dependence becomes less 
drastic. In particular, for U = 2.5t and U = 2t we find a d-wave superconducting region, 
which is, however, smaller than in the frequency- independent calculation. 

In a next step we improve the frequency parameterization of boson propagators by dis- 
cretizing the transfer frequency dependence, rather than assuming a Lorentz form. The 
corresponding stopping scale is plotted in Fig. |3] (line iii). For t' < 0.3t and t' > 0A2t 
the flow of the most singular couplings obtained in the simple Lorentz parameterization 
scheme is essentially reproduced. This is remarkable because, as has already been pointed 
out above, simple Lorentz distributions describe the dependence on transfer frequency well 
only in the small frequency regime. The full discretized frequency parameterization gives a 
slower frequency decay for intermediate to large frequencies for the antiferromagnetic bo- 
son propagator Mii(/o,7r). Furthermore, Mii(/o,0) decays a little faster, see Fig |H This, 
however, does not entail a significantly larger rf-wave superconducting coupling. In the an- 
tiferromagnetic and in the ferromagnetic regime the stopping scale is reduced only by a 
negligible amount. 

While this seems to encourage a simpler Lorentz curve setup over a more involved fre- 
quency parameterization, the former does not allow sign changes of boson propagators in 
their frequency dependence. This restriction is not present in the discretized parameteriza- 
tion. In fact, in the density channel at zero transfer momentum we find a sign change of 
-^ii(^0;0) at a non-zero transfer frequency. The channel decomposition is set up such that 
the boson propagators are positive for transfer frequency zero. Positive values of Kn in the 
s-wave channel correspond to an attractive interaction of density pairs, the flow of which is 
dominated by the local Coulomb repulsion U and can therefore not diverge. Once Ku^Iq, 0) 
becomes negative, i.e. repulsive, it will be enhanced by U. The sign change of the right hand 
side in Eq. (!3T|) is due to the frequency dependence of boson propagators that enter vertex 
and box diagrams via af, which decay in the loop frequency. The biggest contribution to 
this effect comes from the magnetic channel. The bare particle-hole bubble vanishes at mo- 
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FIG. 4: (Color online) Comparison of effective widths A{B) of boson propagators Mii(/o,0), 
Mii(/o,7r), and maxp~^ Mii(/o, p) (from left to right) obtained from Lorentz (black dashed line) 
and discrete (blue line) parameterizations at t' = 0.3t and U = 3t. Here A(i?) is defined by 
B{A,p) = ^B{0,p). Generally, the effective width decreases with Q leading to narrow peaks at 
low RG scales. 

mentum 1 = for frequencies Iq ^ 0. Although this is changed by the f2-regularization we 
expect only minor contributions from the direct graph. The direct graph vanishes for the 
T-scheme, see below. 

We find that a negative value of -ft'ii(/o, 1 = 0) becomes the leading coupling in the 
parameter range 0.3 ^t'/t< 0.42 for U = St. Here, the scattering propagator Ku^Iq, 1 = 0) 
slowly grows to positive values at /q = during the flow, whereas in a certain frequency 
range away from /q = it quickly runs to negative values, compare Fig. [5l The frequency 
Iq'^^{Q) of minimal Kh^Iq, 1 = 0) is roughly proportional to Q, see Fig. [61 We have verified 
that this scattering divergence is robust against refinement of the discretization of exchange 
propagators in frequency and momentum space. 

By decreasing U we again find a growing region where rf-wave superconductivity is the 
leading instability, see the right plot in Fig. |3l The i^-scattering remains leading in a region 
between (i-wave superconductivity and ferromagnetism. 

Up to now we have assumed frequency independent boson-fermion vertex functions, which 
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FIG. 5: Scale dependence of the function Iq i— t- i^ii(Zo,l = 0) for U = 2.5t, t'/t = 0.33. The 
stopping scale is about $7* ~ 0.012t, the initial scale is il.Q = 40i. The function is shown in the 
course of the flow for n/Q^ ^ 1600,400,150,40,20 (left) and w 20,4,2.5,1.5,1 (right) with 

progressive forming of the minimum at non-zero frequency. 




FIG. 6: (Color online) The frequency Zq"^ of minimal i^rii(/o,0) in dependence of the RG scale fl 
for t' = 0.3t and U = 3t (blue line). Due to the finite frequency grid the dependence is a step 
function. For comparison the black dashed line denotes Zq"^ = fi. 

were simply given by scale independent form factors. In a next step we evaluate the full flow 
equations given in Sec. llVt where in addition to boson propagators we also discretize the 
frequency dependence of 7i 2"*, 7^, and 7^. Their frequency dependence is plotted in Fig. [71 
The s-wave channels deviate only by a small amount from unity (which is the normalization 
at frequency zero). For large frequencies, attractive channels saturate at a value below 1 
and repulsive channels above 1. In the d-wave superconducting channel the boson- fermion 
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FIG. 7: (Color online) Frequency dependence of boson-fermion vertex functions at 0* ~ 0.037t for 
U = 3t and t' = 0.3 at Van Hove filling: (a) 75"^'** , (b) , (c) 7^ , and (d) 72*^'*^. The blue curves 
correspond to a = and the green dashed lines to a = vr. 



vertex function decays to zero. Despite the frequency dependence of the boson-fermion 
vertex functions we find only minor changes in the stopping scale, see symbols in Fig. [3l 
Most importantly, the i^-scattering process remains leading in the same parameter region 
as before. 

We cannot exclude the possibility that the /T-scattering singularity with sign change in 
the density channel is an artifact of the channel decomposition that is not present in the full 
level-2 truncation of the RG. However, our result is stable against the inclusion of frequency 
dependent fermion-boson vertex functions. For determining singularities in the RG flow 
the dependence of the vertex function on non-transfer frequencies seems negligible. This 
dependence may become important in a more quantitative analysis of observables. 



Even if Kil^"^^, 0) is not the leading coupling at the stopping scale, its sizeab 
to the effective interaction cannot be disregarded. In the conventions of Ref. 



e contribution 



13 it reads 



Vk[*] = -\j^l K{l)Y{l)Y{-l) (48) 

with density operators Y{1) = f dp \E'o-(p)^'o-(p + 0- Suppose that Ji (/q, 0) has a strong 
negative peak at /q"*^. In the following we neglect contributions from other channels and 



28 



assume a mean-field approximation 

i^(0 = -4A5(l)5(|/o|-/r) (49) 

at the critical scale with A > 0. Then the i^-scattering singularity corresponds to a repulsion 
of time modulated densities 

V^^[^]=A J2 n{x,T)n{x',T')e<"^'~^'^ (50) 

x,x',t,t' 

where x denotes lattice sites and r is the imaginary time. For /q*"^ = 0, mean field theory 
produces only a shift of the chemical potential. In contrast, /q*"^ > 0, which is the situation 
that we encounter, yields off-diagonal frequency terms in the two-point function and hence 
changes its structure in a non-trivial way. A detailed analysis of this scattering singularity, 
in particular regarding its influence on an RG flow below the stopping scale is left for future 
work. 

B. Influence of Imaginary Exchange Propagators 

Unlike in RG setups where frequency dependences are neglected, the frequency depen- 
dent interaction vertex acquires an imaginary part during the flow. By symmetry, in our 
parameterization particle- hole exchange propagators Mn, Ku are real- valued whereas Dn, 
D22 in the particle-particle channel are not. Still, lmDmm{io = 0,1) = is imposed by 
symmetry. 

Perturbation theory and RPA resummation yield a singular behavior in the first frequency 
derivative of ImDu, originating from an asymmetry in the density of states about zero 
energy, see Eq. (143|) . We have examined the influence of this singularity in the RG calculation 
at the data point t'/t = 0.3, U/t = 3. Here the stopping scale is relatively low such that an 
effect, if present, should be visible. Furthermore, a strong scattering singularity has been 
encountered in this parameter region and we study how it is influenced. 

The fact that ImDmmiO exactly vanishes at /q = already suggests that ImD^m pos- 
sibly has limited impact. This is confirmed by carrying out the RG calculation: Although 
{dig lmDmm){lo = 0, 1) shows the expected singular behavior, ImDmm remains small during 
flow. At the stopping scale, |ImDii(/)| < 5t with the maximum at 1 = tt and Iq ^ 10^2^,, as 
well |Im 1)22(0 1 < 0.05t. 
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FIG. 8: (Color online) Influence oiluiDmm. on the flow of the most singular couplings for frequency 
dependent interaction vertex. The flow with slightly lower stopping scale takes into account Im Dn 
and lmD22, the second one neglects them. Line conventions as in Fig. [3j 

As a consequence, these quantities have little influence on the flow of the leading cou- 
plings, see Fig. |H] where the flows including and neglecting the imaginary parts of the SC 
channel are compared. In particular, the fC-scattering singularity remains unchanged. Al- 
though we have not performed a detailed scan of the whole parameter space we expect 
ImDmm to be of minor influence. 

C. Comparison to T-Scheme 

We need to check whether the /^-scattering singularity is an artifact of the Q- 
regularization. To this end, we repeat the calculation in the temperature flow scheme^, 
which uses temperature T as the scale parameter. This scheme is used frequently in RG 
studies of the Hubbard modeU^i^, so far mostly in the static approximation. For our purpose 
the frequency dependence of the interaction vertex needs to be taken into account as well. In 
terms of rescaled fields, the free propagator in the T scheme reads Ct{p) = T^^'^[ipo — e{p)]~^ 
and fermionic Matsubara frequencies Pq G 7rT(2Z — 1) are discrete. 

The flow equation ([3]) can be rewritten by setting up a channel decomposition as before. 
Exchange propagators now carry discrete bosonic frequencies Iq € ttT TL. We project non- 
transfer frequencies to the lowest possible frequency value, i.e. po = ±vrT, and symmetrize 
over the sign. Self-energy effects and imaginary parts of the interaction vertex are neglected. 
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FIG. 9: (Color online) Temperature flow at Van Hove filling. Left: Flow of most dominant insta- 
bilities. The stopping temperature here is T* ~ 0.22t. Line conventions as in Fig. [3l Right: Flow 
of Kii{lo,0) for the lowest possible frequency values Iq/ttT = 0,6,4,2 (from top to bottom). The 
scattering singularity shows at the lowest non-zero transfer frequency. 

We calculate the flow for 0.3 < t' /t < 0.35, i.e. in the parameter region where a strong 
scattering process has been found in the f2-scheme, and for several Hubbard parameters U. 
The initial condition of the flow is determined from second order perturbation theory at a 
high temperature Tq = lOOt. Technical details of the calculation are given in Ref. 
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Firstly, we observe that the T-flow diverges faster than the f2-flow in the sense that, 
when comparing scales as energy variables, T* is much larger than (whereas in the static 
approximation ^ Q^). We interpret this as a consequence of the specific form of transfer 
frequency dependence, which is quite different in both schemes. Whereas we observe that 
exchange propagators are monotonic functions of the scale parameter in the T-scheme, we 
find non-monotonic behavior for the f2-scheme. Here, exchange propagators at a specific 
transfer frequency Iq during the flow usually first grow for Q > Iq and then shrink for 
Q < Iq. Consequently, as a function of frequency, exchange propagators in the f2-scheme are 
typically broader than in the T-scheme and hence closer to the approximation of frequency 
independent exchange propagators. 

Furthermore, the zero-momentum scattering exchange propagator Kii{Iq, 0) again shows 
two different behaviors: At zero frequency it slowly flows towards positive values while for 
the smallest non-zero exchange frequency, i.e. Iq = ±27rT, it strongly grows to negative 
values, compare Fig. [91 This corresponds to /q"*^ = 27rT in the T-scheme and is consistent 
with the approximately linear decrease of /q"*^ as a function of scale in the f2-scheme. 
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We have not found parameter values where this scattering process becomes the dominant 
couphng in the T-flow. However, it is present and gets strong in this scheme as well. For 
that reason we do not consider it an artefact of a particular regularisation. 

D. Consideration of Imaginary Self-energy 

A non-trivial flow equation for the frequency dependent self-energy can be constructed 
already with a static approximation of the interaction vertex—. Knowledge about the full 
frequency dependent interaction vertex now allows a detailed study of the frequency de- 
pendent self-energy. Analysis of in second order perturbation theory^l for a system 
at van Hove filling and with square Fermi surface shows a logarithmic singularity in the 
first frequency derivative at momentum p = (0,7r). It turns out that this singularity gets 
enhanced in the RG setup and results in a suppression of the interaction vertex during the 
flow. This suppression can become rather strong. 

We calculate the flow of the imaginary part of the self-energy, for which most singular 
behavior is expected, 

lmt{p) = - j d/ \mS{l + p)(^ReD + h<+h'I)i^{l) + ReD22{l)fl{\ + v)^. (51) 

Imaginary parts of exchange propagators in the pairing channel are neglected here since 
they remain small^^. Exchange propagators are discretized in frequency and momentum 
space as before. We then trace the function uj ^ ImS(u;,p) by frequency discretization 
for several momenta p. Note that a simple parametrization of Im S using only a linear 
(momentum-dependent) frequency term is insufficient and leads to an artificial suppression 
of the fiow. 

It turns out that the presence of the full propagator (i.e. including Im S) instead of the 
bare one makes the evaluation of the fiow equations relatively time-consuming since all 
loop frequency-momentum integrations then need to be performed numerically at each RG 
step. At low scales the involved triple integration in the full fiow equations was not done 
because of long computing times. As an approximation, for parameter values with fiow to 
low scales, we restrict to a momentum independent feed-back of ImS(p) to the fiow. This 
approximation again allows factorization of the propagator and largely simplifies numerics 
by permitting computation of loop momentum integrals beforehand, similar to the strategy 
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FIG. 10: (Color online) Comparison of stopping scales for different RG setups at Van Hove filling. 
Top-down: (i) frequency dependent vertex disregarding self-energy, (ii) frequency independent 
vertex disregarding self-energy, (iii) frequency dependent vertex including Im S(a;, (0, vr)). Line 
conventions as in Fig. [31 
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FIG. 11: (Color online) Dependence of the Z-factor scale derivative on different frequency pa- 
rameterizations of exchange propagators. Compared to the correct scale derivative (computed 
from frequency discretized exchange propagators, -|- marks), the parameterization by a sum of two 
Lorentz curves (x marks) yields precise results. A single Lorentz curve determined from the small 
frequency behavior (circles) overestimates the flow. 

of Sec. IVI C[ Thus, we replace ImS(a;,p) — )■ ImS(a;, (0,7r)) in the right hand side of the 
flow equations. We have chosen the momentum (0, tt) for the self-energy feed-in since the 
saddle point region, for systems at Van Hove filling, mainly drives the flow and should receive 
special attention. Because in this region ImS also shows its most singular behavior, this 
approximation overestimates the suppression of the flow. 
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The resulting stopping scale is shown in Fig. [TO] for initial interaction U = 3t. We observe 
three main differences as compared to the frequency dependent setup disregarding self-energy 
effects. First, the stopping scale is significantly lowered and drops drastically in the 
parameter region of competing pairing and ferromagnetic ordering tendencies. Furthermore, 
a region of dominant d-wave pairing is recovered. On the contrary, the i^-scattering is 
weakened and no longer becomes dominant. 

The flow of the self-energy is very sensitive to proper frequency parameterization of 
exchange propagators. In the calculation of the flow as described above we additionally de- 
termine the scale derivative of the 2'(o,7r) factor with alternative frequency parameterizations 
of exchange propagators and compare to the correct result in Fig. [TTJ This test does not 
take into account an error propagation because at each scale the true self-energy is used to 
calculate scale derivatives. We observe that the parameterization of exchange propagators 
by a sum of two Lorentz curves produces quantitatively accurate results. 

A more detailed study about self-energy flows as well as a discussion about the neglect 



of the momentum dependence of Im S at low scales can be found in Ref. 
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VIII. CONCLUSION 



We have presented a functional RG calculation for the two-dimensional Hubbard model 
at Van Hove filling. All calculations use a channel decomposition for the interaction vertex, 
which provides an efficient parameterization of momenta and frequency. We studied the 
effective two-fermion interaction and compared quantitatively the effect of different approx- 
imations of the frequency dependence. 

In general, including frequency dependences of vertices in an RG flow with bare prop- 
agators raises the (pseudo-) critical scale compared to frequency independent calculations. 
The higher scale is accompanied by a reduction of the d-wave superconductivity coupling 
that becomes leading in a smaller region of the {U, t')-parameter space. In particular, for 
U = 3t we did not find a d-wscve instability at Van Hove filling. Lower initial repulsive cou- 
plings U re-establish a rf-wave instability and the approximation of neglecting the frequency 
dependence leads to less drastic effects. 

In a second calculation we have included the imaginary part of the fully frequency dis- 
cretized self-energy, evaluated at the Van Hove points. The divergence of the corresponding 
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Z factor then leads to a suppression of the flow of the frequency dependent interaction 
vertex. The stopping scale is reduced and can drop to rather small values, notice that the 
evaluation of the self-energy at the Van Hove points overestimates this reduction. Besides 
the change in the stopping scale, the flow qualitatively agrees well with previous flows that 
neglect frequency dependences at all. Again we flnd regions of dominating antiferromag- 
netism, (i-wave superconductivity and ferromagnetism, whereas couplings in the scattering 
channel remain subleading. In this view, neglect of frequency dependences and self-energy 
contributions at the same time seems to be a rather good approximation for the full RG flow. 
The reason for this is not fully understood yet. It is, however, clear from the lowest equation 
in the RG hierarchy that keeping a frequency-independent self-energy in a calculation with 
a frequency-dependent vertex cannot be exactly correct. In fact, such an approximation 
violates the Ward identity between the F- vertex and S. Including the frequency dependence 
in our case restores the results of the static approximation qualitatively. 

For future applications of the functional RG it is important to develop feasible approxi- 
mations of the frequency dependence. The simplest parameterization with a single Lorentz 
curve for each boson propagator already gives reasonable results for the flow of dominant 
couplings in the vertex function. Since this ansatz does not correctly describe intermediate 
to large frequencies, which contribute more to the flow equation for the self-energy than in 
the vertex function, the error for the Z-factor becomes large. With two Lorentz curves each 
boson propagator can already be described well enough to compute the imaginary self-energy 
and the leading instabilities of the interaction vertex reliably. 

The benchmark for testing frequency parameterizations was set by a full discretization 
of frequency and momenta for each boson propagator. Imaginary contributions to the in- 
teraction vertex seem to have little impact on the RG flow. While single Lorentz curves 
do not allow sign changes in the frequency dependence, we flnd indeed such a sign change 
in the density forward scattering channel in the discretized flow. The corresponding cou- 
pling at transfer momentum zero and transfer frequency /q"^ can in principle diverge and 
becomes large in a region between superconductivity and ferromagnetism. The frequency 
is roughly proportional to the RG scale. 

While in the f2-regularization scheme the coupling of this this flnite frequency K- 
scattering becomes leading in a flow with bare propagators, this is not the case in the 
temperature regularization scheme or when the imaginary self-energy is included. In any 
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case we find a large contribution to the effective interaction at finite scales. Its influence on 
low scales and on observables is left for future work. 

All calculations were performed using a channel decomposition of the level-2 truncation, 
which made the parameterization of momenta and frequencies more manageable. We gen- 
eralized the method to be able to describe non-transfer frequencies. Unlike non-transfer 
momenta they are not expanded in orthonormal functions. Instead we derived flow equa- 
tions for fermion-boson vertex functions, that is, frequency-dependent Yukawa couplings. It 
turned out that their frequency dependence does not significantly change results for the RG 
flow of dominant couplings in the Hubbard model. 

We thank A. Eberlein, J. Ortloff, and W. Metzner for fruitful discussions. This work was 
supported by the DFG research unit FOR 723. 
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